A new family of time-space harmonic polynomials with respect to 

Levy processes 



^ ■ E. Di Nardo * I. Oliva t 

O 

(N : April 27, 2012 

u 

< 

^ ■ Abstract 

(N 

By means of a symbolic method, a new family of time-space harmonic polynomials with 
respect to Levy processes is given. The coefficients of these polynomials involve a formal 
expression of Levy processes by which many identities are stated. We show that this fam- 
ily includes classical families of polynomials such as Hermite polynomials. Poisson-Charlier 
polynomials result to be a linear combinations of these new polynomials, when they have the 

■ property to be time-space harmonic with respect to the compensated Poisson process. The 
more general class of Levy-Sheffer polynomials is recovered as a linear combination of these 
new polynomials, when they are time-space harmonic with respect to Levy processes of very 

_ general form. We show the role played by cumulants of Levy processes so that connections 

^ ■ with boolean and free cumulants are also stated. 
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O I 1 Introduction 

A family of polynomials {P{x, t)}t>o is said to be time-space harmonic with respect to a stochastic 
process {Xt}t>o if E[P{Xt,t) \ ^s] = PiXs,s), for all s < t, where '^s = cr {Xr ■ t < s) is the 
^ ■ natural filtration associated with {Xt}t>o- For random walks {Xn}n>o, Neveu [lOj characterizes 

^ ■ the family of time-space harmonic polynomials as the coefficients of the Taylor expansion 



£[exp{=X„)] 'kl 

in some neighborhood of the origin. If {Xn}n>o is replaced by a Levy process {Xt}t>o, the left- 
hand side of (II. ip is the so-called Wald's exponential martingale [9]. The usefulness of time-space 
harmonic polynomials with respect to Levy processes is that the stochastic process {P{Xt,t)} is 
a martingale, whereas {Xt} does not necessarily have this property. 

The Wald's exponential martingale is well defined only when the process admits moment 
generating function i?[exp{2;Xt}] in a suitable neighborhood of the origin. Different authors have 
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tried to overcome this gap by using other tools. Sengupta [15j uses a discretization procedure 
to extend the results proved by Goswami and Sengupta in [8]. Sole and Utzet [T7| use Ito's 
formula showing that time-space harmonic polynomials with respect to Levy processes are linked 
to exponential complete Bell polynomials yj. The Wald's exponential martingale (jl.ip has been 
recently reconsidered also in [16], but without this giving rise to a closed expression for these 
polynomials. 

In this paper, by using a symbolic method, known in the literature as the classical umbral 
calculus, we give a new family of time-space harmonic polynomials, which could be easily imple- 
mented in any symbolic software, see [5] as example. Thanks to the results in [1], we show that 
this new family includes and generalizes the exponential complete Bell polynomials. 

The classical umbral calculus we use consists essentially in a moment symbolic calculus, since 
its basic device is to represent an unital sequence of numbers by a symbol q, named umbra, i.e. to 
associate the sequence 1, ai, a2, . . . to the sequence 1, a, a^, . . . of powers of a through an operator 
E that looks like the expectation of random variables The n-th element of the sequence is 
the n-th moment of a. As a matter of fact, an umbra looks the framework of a random variable 
with no reference to any probability space. 

In this paper, we define an operator S[-|a] that acts like the well-known conditional expectation 
of random variables. The umbral version of Levy processes we propose takes into account their 
infinite divisible property. As corollaries, many identities are given on the coefficients of this family 
of time-space harmonic polynomials with respect to random walks and Levy processes. Moreover, 
this expression allows us also to emphasize the role played by cumulants and to include boolean 
and free cumulants [3l[6]. As example, we show that this new family includes Hermite polynomials 
which are time-space harmonic with respect to Brownian motion. We prove that Poisson-Charlier 
polynomials are linear combinations of the introduced polynomials when they are time-space 
harmonic with respect to compensated Poisson processes. We show that the more general class 
of Levy-Sheffer polynomials [H] are linear combinations of the introduced polynomials and 
we characterize a general form for the associated Levy processes. 

The paper is structured as follows. Section 2 is provided for readers unaware of the classical 
umbral calculus. Let us underline that the theory of the classical umbral calculus has now reached 
a more advanced level compared to the elements here resumed. We have chosen to recall termi- 
nology, notation and the basic definitions strictly necessary to deal with the object of this paper. 
In Section 3 the new notion of conditional evaluation with respect to umbrae is introduced, which 
is the key to characterize time-space harmonic polynomials in terms of umbrae. Examples and 
applications are introduced in Section 4, where special emphasis is devoted to the role played by 
cumulants in the expression of these polynomials. 

2 The classical umbral calculus 

In the following, terminology, notation and some basic definitions of the classical umbral calculus 
are recalled. We skip any proof: the reader interested in deeper analysis is referred to the papers 

[21 S]. 

Let M[x] be the ring of polynomials with real coefficients in the indeterminate x. The classical 
umbral calculus is a syntax with an alphabet A = {a,/3,7, . . .} of elements, called umbrae, and a 
linear functional E : M.[x][A\ — > M.[x], called evaluation, such that £'[1] = 1 and 

E[x'' a' ■■■ 7*=] = E[a'] E[p^] • • • E[-/''] (uncorrelation property) 
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where a, /3, . . . , 7 are distinct umbrae and n,i,j, . . . ,k are nonnegative integers. 

A sequence {a„}„>o G with oq = 1, is umhrally represented by an umbra a \i E[a^] = an, 
for aU nonnegative integers n. Recall that a„ is called the n-th moment of a. An umbra is scalar 
if its moments are elements of M while it is polynomial if its moments are polynomials of 
Special scalar umbrae are: 

i) the augmentation umbra e, such that E[e^] = 5o,nQ for all nonnegative integers n; 

ii) the unity umbra n, such that E[u'^] = 1, for all nonnegative integers n; 
Hi) the i?eZZ umbra /3, whose moments are the Bell numbers; 

iv) the singleton umbra x such that E[x\ = 1 and E[x^] = 0, for all integers n > 1. 

The core of this moment symbolic calculus consists in the definition of the dot-product of two 
umbrae, which is fundamental both in the construction of time-space harmonic polynomials and 
in their applications. We recall in short the steps necessary to give this definition. 

First let us remark that in the alphabet A two (or more) distinct umbrae may represent the 
same sequence of moments. More formally, two umbrae a and 7 are said to be similar when 
^[a"] = £^[7"] for all nonnegative integers n, in symbols = 7. Therefore, given a sequence 
{on}n>0; there are infinitely many distinct, and thus similar umbrae, representing this sequence. 

Now, define the symbol n.a representing a' + a" + • • • + a'", where {a', a", . . . , a'"} is a set 
of n uncorrelated umbrae similar to a. The symbol n.a is an example of auxiliary umbra. In a 
saturated umbral calculus, the auxiliary umbrae are treated as they were elements of A [H]. The 
umbra n.a is called the dot-product of the integer n and the umbra a. Its moments are [2]: 

i 

qi{n) = E\{n.ay] = ^(n)fcSj_fc(ai, 02, . . . , ai^k+i), (2.1) 
k=l 

where (n)^ is the lower factorial and Bi^k are the exponential partial Bell polynomials [Tj. 

In (|2.ip . the polynomial qi{n) is of degree i in n. If the integer n is replaced by t S M, in (|2.ip 
we have qi{t) = X]fc=i(^)fe-^i,fc(^i' (^2-, ■ ■ ■ ■, Oi-fe+i)- We denote by t.a the auxiliary umbra such that 
E[[t.ay] = qi{t), for all nonnegative integers i. The umbra t.a is the dot-product of t and a. 
Among its properties, we just recall the distributive property: 

(t + s).a = t.a + s.a', s,t£R (2.2) 

where a' = a. In particular in ()2.ip we can replace n with —t by obtaining the auxiliary umbra 
—t.a with the remarkable property 

-t.a + t.a' = e, (2.3) 

where a' = a. Due to property (|2.3p . the umbra —t.a is named the inverse umbra of t.a. H 

Let us consider again the polynomial qi{t) and suppose to replace t by an umbra 7. The 
polynomial qi^j) is an umbral polynomial in M[x][^], with support supp {qi{^)) = {7}. Recall that 
the support supp (p) of an umbral polynomial p G M[x][^] is the set of all umbrae occurring in 
it. The dot-product of 7 and a is the auxiliary umbra 7.0 such that £^[(7.0)*] = E[qi{^j)] for all 

^The symbol Sij denotes the Kronecker delta, that is 5ij = 1 if i — j, otherwise Sij = 0. 

^Since —t.a and t.a are two distinct symbols, they can be considered uncorrelated, therefore —t.a + t.a' = 
—t.a + t.a = e. When no confusion occurs, we will use this last similarity instead of (|2.3p . 
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nonnegative integers i. Two umbral polynomials p and q are said to be umbrally equivalent if 
E[p] = E[q], in symbols p q. Therefore equation (j2.ip . with n replaced by an umbra 7, can be 
written as qi{'^) ^ (t-Q;)' ^ Sfc=i(7)A:A,fc(oii • • • j Q^i-fc+i)- Special dot-product umbrae are the 
a-cumulant umbra x-ct and a-partition umbra (3. a, that we will use later on. Recall that if 0,7 
and ?7 are uncorrelated umbrae, then 

(a + r7).7 = a.7 + r/.7. (2.4) 

3 Time-space harmonic polynomials 

Denote by X the set X = {a}. 

Definition 3.1. The linear operator E{-\a) : — > W^X] such that 
i) E{1 I a) = 1; 

a) E{x"^a^^^5^ • • • I a) = x'^a"'E[y'^]E[5^] ■ ■ ■ for uncorrelated umbrae 0,7, 5, . . . and for non- 
negative integers m, n, i, j, . . . 

is called conditional evaluation with respect to a. 

In other words, Definition 13 . 1 I savs that the conditional evaluation with respect to a handles the 
umbra a as it was an indeterminate. The proofs of the following Propositions are straightforward 
taking into account Definition 13. 1[ 

Proposition 3.2. If a ^ A and p € M[j;][^] with a supp(p), then E{p \ a) = E\p\. 
Corollary 3.3. If a e A and p e W\x\[A], then E[E{p \ a)] = E[p\. 

This last corollary brings to light the parallelism between the conditional evaluation E{-\a) 
and the well-known conditional expectation in probability theory [7] . As it happens in probability 
theory, the conditional evaluation is an element of M[a;][^] and, if we take the overall evaluation 
of E{p I a), this gives E[p\. 

The conditional evaluation with respect to the auxiliary umbra n.a is such that E[{n + 
\).a I n.a] = E{n.a + a' \ n.a) = n.a + E[a'], with a' an umbra similar to a. By similar ar- 
guments, for all nonnegative integers n and m we have 

E{\{n + m).a]^ \ n.a) = E{[n.a + m.a']^ \ n.a) = ^ ( ^ j {n.ay E[{m.a')'^^'^], 

and by taking the evaluation of both sides we recover E{[{n + m).a]'^) due to distributive property 
(j2.2p . Therefore, for f > we define the conditional evaluation of t.a with respect to the auxiliary 
umbra s.a, with < s < i such as 

E[{t.a)'' I s.a] = (^) {s.ayE{[{t - s).af-^). (3.1) 

Definition 3.4. Let {P{x,t)} £ M[x] be a family of polynomials indexed by t >0. P(x,t) is said 
to be a time-space harmonic polynomial with respect to the family of auxiliary umbrae {qit)}t>o 
if and only if E [P{q{t),t) \ q{s)] = P{q{s), s) for all < s < t. 
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Theorem 3.5. For all nonnegative integers k, the family of polynomials 

Qk{x, t) = E[{x - t.a)^] E M[x] (3.2) 
is time-space harmoni^ with respect to {t.a}t>o. 

Proof. From (|3.2p . by applying the Unearity property of the evaluation E, we have 

k 



Qkix,t) = j2(^)^''-'E 



[{-t.ay] (3.3) 



j=0 

for all nonnegative integers k. Thanks to ()3.ip and ()2.2p . we have 

k 



E {Qk{t.a, t) I s.a) = ^ E[{t.af-^ \ s.a] E[{-t.ay] 



j=0 ^■'^ l.j=o 
By suitably rearranging the terms, we have 

/ , \ ( k-j 



E{Qk{t.a,t)\s.a) = {s.ay I^J^I^^ i ^ 



Corollary 3.6. IfQk{x,t) = ELo (*) th 



□ 



en 



i) q^l'\t) = {^^E[{-t.af-^, forty Q and j = 0,1, . . . ,k; in particular qf\t) = 1; 

ii) qf\{)) = /or J = 0, 1, . . . , - 1 and q[^\o) = 1. 

In particular we have Qk{x,0) = x^ for all nonnegative integers k. 

Proof. Property i) follows from (|3.3p . Property zzj follows by observing that (/^'^^(O) = (J)£'[(0.q)'^~-'] 
^[gfc-i] for j = 0, 1, . . . , A: - 1 and ^(0) = ^(e^) = 1. □ 

The sequence of polynomials {Qk{x,t)} is umbrally represented by the polynomial umbra 
X — t.a. We call x — t.a the time-space harmonic polynomial umbra with respect to t.a. 

Proposition 3.7. The time-space harmonic polynomial umbra x — t.a is the Appell umbra of 
—t.a. 



'When no confusion occurs, we will use the notation x—t.a to denote the polynomial umbra —t.a+x = x+{—t.a). 



5 



The result follows by the definition of Appell umbra given in [3]. In particular Proposition 13.71 
means that the sequence of polynomials {Qk{x,t)} is an Appell sequence, that is 

—Qkix,t) = kQk-i(x,t), for all integers A; > 1. 
dx 

Let us observe that every linear combination of polynomials {Qkix,t)}k>i is a time-space 
harmonic polynomial with respect to {t.a}i>o. 

Remark 3.8. Note that the family of auxiliary umbrae {t.a}t£i, with / C M+, is the umbral 
counterpart of a stochastic process {Xt}t^j such that E[X^] = E[{t.a)'^] for all nonnegative 
integers k. This stochastic process is a Levy process, see next section and in particular Remark 
14.51 for a parallelism between the definition of t.a and the infinite divisible property of a Levy 
process. Therefore, the polynomials {Qk{x,t)} are time-space harmonic with respect to Levy 
processes. If the moments of {Xt}t£i are defined only up to some finite m, the representation 
()3.2p still holds up to m, because it involves moments of order less or equal to m, see also next 
remark. 

Remark 3.9. The generating function of an umbra a is the formal power series f{a,z) = 
X^n>o'^"7n" ^ ^[^lli-^]] whose coefficients are the moments of the umbra, see [2] for more details. 
Formal power series allow us to work with generating functions which do not have a positive radius 
of convergence or having undefined coefficients The generating function of the time-space 

harmonic polynomial umbra x — t.a is 

fix - t.a, z) = ^-^^ = Qkix, tf^. (3.4) 

By replacing x with t.a in (13. 4p . we recover the Wald's exponential martingale (II. ip . Equality of 
two formal power series is interpreted as the equality of their coefficients, so that E[Ri:{Xt,t)] = 
E[Qk{t.a,t)]. 

Also the Wald's identity Y.k>o ^l^k{Xt,t)]z'' /k\ = 1 follows from ([TI]). Therefore, the 
sequence {£'[i?fc(Xt, t)]}fc>o G M is umbrally represented by the augmentation umbra e, with 
/(e, z) = 1, for all t > 0. But this is exactly what it happens when in the polynomial umbra 
X — t.a we replace x with t.a. 

The following corollary specifies the dependence of the coefficients of Qk{x,t) in ()3.2p on the 
umbra a. 

Corollary 3.10. If {an} is the sequence umbrally represented by the umbra a and {Qki^^t)} '""e 
time-space harmonic polynomials with respect to {t.a}t>o, then Qk{x,t) = '^ji^Qcf'jt^x^, with 

'k 

d 



) ^ d,(-l)2'W+^s[/(A),i]a?a^^--- 

where s[l{X),i] are the Stirling numbers of first kind, the sum is over all partitionEx = ,2""^ , . . .) 
of the integer k-j and dx = H/irM ■ ■ ■ (l!)''i (2!)''2 . . . )_ 



^Recall that a partition of an integer i is a sequence A — (Ai, A2, . . . , Am), where Aj are weakly decreasing positive 
integers such that X^JLi ~ The integers \j are named parts of A. The length of A is the number of its parts 
and will be indicated by /(A). A different notation is A = (V^ , 2^^ ,■■■), where rj is the number of parts of A equal 
to j and ri + r2 + ■ ■ • = Note that rj is said to be the multiplicity of j. We use the classical notation A h i to 
denote "A is a partition of i" . 
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Proof. Equation ()2.ip . with n replaced by t, can be restated as E[{—t.ay] = ^xhi ^xi't)i(x)^V^2 
see [3]. In particular we have ^[(-t.a)'=-J] = EAhfc-j cIa s[l{X),i] {-ly+^^^^U^ a['al^- 

where s[/(A),i] is the i-th Stirling number of the first kind. From Corollary 13.61 we have 



j=0 \xhk-j 



2 ■■■Y,i-'^y~'''^'^4m^¥\x' (3.5) 



i=0 



and the result follows by suitably rearranging the terms in (|3.5p and by observing that s[l{\),i] = 
for i > /(A). □ 

In the following, assume Qk{x,t) in (j3.2p such that Qk{x,t) = X] j=o ^j'^'' (^) ^^^"^ denote by 
{ttn} the sequence of moments umbrally represented by the umbra a in (|3.2p . 

Proposition 3.11. We have qf\t - 1) = Y.Lj (j^f ^(0 <^i-j- 



Proof. Prom Corollary [Ml we have q) '{t - 1) ~ - ^ {'^){-t.a + 

(J) Es=o C7){-t.a)^-^-'a' ~ Q) «*"^ - The result follows by taking the eval- 

uation E of both sides. □ 

Proposition 3.12. We have at = q^Q\t - 1) - Y^'^j^o ifH^) f^j- 

Proof. By using Proposition 13. IH we have qQ''\t — 1) = akq^j^\t) + X]j=o We have 
<£\t - 1) - E -=0 ^^qf\t) = akqi'Ht) + EjZo mf\t) - E -=0 mf^it) = cikqi'\t). The result 
follows by observing that q^^^ (t) = 1. □ 



Theorem 3.13. A polynomial P{x,t) = Ej=o^'i(^) degree k for all t > 0, is a time-space 



harmonic polynomial with respect to {t.a}t>o if and only if 

k 

M*) = E ('•) ^'^(O) E[{-t-c^n /or J = 0, ... , k. (3.6) 
i=j ^-^^ 

Proof. Assume P{x,t) = E j=o ^'i (*) ^ polynomial whose coefficients satisfy (13. 6p . Then we 
have 

E I E (')l'.(0)i?[(-t.«)^-^] i = Ef^WE ^^^[(-i-«)n = Y.Pi{0)E[{x-t.an 
j=0 [i=j ^-^^ J j=0 i=0 i=o 

(3.7) 

As P(x,t) is a linear combination of {Qfc(a;,t)}, then P(x,t) is a time-space harmonic polynomial 
with respect to {t.a}t>o- Viceversa if P{x, t) = Yl^=oPj(^) ^ time-space harmonic polynomial 
with respect to {t.a}t>o, then P{x,t) = Yli=o'^i'^[ix — t.aY], with {q} S M. Therefore, from (j3.7p 
and for j = 0, . . . , A; we have Pj(t) = Ei=j t.a)*"-']. So (|3.6p follows by observing that, 

when t is replaced by 0, we have pj{0) = j^Lj Ci E[{-0.ay-^] = J2i=j Q) d E[e'-^] =Cj. □ 

Corollary 3.14. If P{x,t) = ^^j^oPjit) x^ is a polynomial of degree k for all t > 0, then there 
exists an umbra a such that P{x,t) is a time-space harmonic polynomial with respect to {t.a}t>o- 
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4 Applications and examples 

4.1 Discrete Case 

When the parameter t is replaced by a nonnegative integer n, the coefficients qj^\n) of time-space 
harmonic polynomials satisfy further properties thanks to the umbral representation (j3.2p . To 
keep the length of the paper within bounds, we just show some of them. 

Proposition 4.1. We have qf\n) + Y!1=j+i Q E?=i = 0- 

Proof. FromCorollary[M]Et,+i T.ti^'~' ^ E?=iEIj+i O'"^' ( ■) (-^•«)'-^ Since 
Et,+i Q)a*-^'(')(-^-«)'"' ^ (?) + "O'--^' - (-^.a)'^--'] =^ qf\l - 1) - 'zf (/), the result 

follows by observing that ^t.+i E?=i »'-'ql'\l) ^ E?=i {ifH^ " 1) " (O) ^ 9^(0) - 
qf\n) ~ -gf V). since qf\o) ~ (J) (-O.a)''-^ ~ 0, for all k ^ j. □ 

Corollary 4.2. W^e /lawe gS''^(n) + Y.?=i Ej=i ^(0 = 0. 

Remark 4.3. Let us observe that the family of umbrae {n.a}n>o corresponds to a discrete 
martingale {Xn}n>o with Xq = and independent and identically distributed difference sequence 
with zero mean. Recall that the difference sequence associated to {Xn}n>o is a sequence of 
random variables {Mn}n>o such that Mq = Xq = and M„ = X„ — for all nonnegative 

integers n. The umbra n.a generalizes X„ = Mi + M2 + • • • + M„. Suppose to remove the identical 
distribution hypothesis on {Mn}n>o '■ in umbral terms, the martingale {Xn}n>o corresponds to 
the umbra ai+a2 + - ■ ■+Oin, where the umbrae {ai , . . . , ««} are not necessarily similar. The time- 
space harmonic polynomials E[{x—n.a)''] need to be replaced by E[{x — l.(ai + 02 + • • • + ctn))^]- 
The properties stated up to now can be recovered by similar arguments. 

4.2 Cumulants 

Any umbra is a partition umbra (cf. [2]). This means that if {a„}„>o is a sequence umbrally 
represented by an umbra a, then there exists a sequence {hn}n>i umbrally represented by an 
umbra k^, such that a = fi.n^. In terms of generating functions we have /(a, z) = exp[/(K^, 2:) — !], 
so that {hn}n>i is the sequence of formal cumulants of {a„}n>cll. The umbra k„ is called a- 
cumulant umbra and we also have = x-a, where x is the singleton umbra. 

Proposition 4.4. For the sequence of polynomials {Qk{x,t)} umbrally represented by the time- 
space harmonic polynomial umbra x — t.a, we have 

Qk{x,t) =Yk{x + hi,h2,...,hk), (4.1) 

with Yk exponential complete Bell polynomials and {hn} the sequence of cumulants of —t.a. 

Proof. We have E[{f3.^)^] = Yk(gi,g2, ■ ■ ■ ,gk) with gn = -£■[7"], for all nonnegative integers n [2]. 
Therefore we will prove (14. if we show that x + .(—t.a) = f3.j for some polynomial umbra 7. 
Choose as umbra 7 the umbra K(^^ u-^+K(^_f a) with K(^x.u) the cumulant umbra of x.u and 

^In the ring of formal power series, given a sequence {a„}„>i, its sequence of formal cumulants {hn}n>i is such 
that 1 + E„>i auz"/n\ = exp (e„>i /inz^/nl) . 
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the cumulant umbra of —t.a. Let us recall that the disjoint sum of two distinct umbrae 61+62 is 
an auxiliary umbra whose n-th moment is -^[^i ] + -^[52 ] for all integers n > 1. We have 

nr/ ; \m ( X + hi n = 1 

^[('^(-«)+'^(--)) ] = I /.„ n>l 

with {/in} the sequence of cumulants of —t.a. The result follows since x + l.a) = f3.H(^x.u) + 

Remark 4.5. As f{f3.Ka,z) = exp[f(Ka,z) — 1], the auxiliary umbra f3.n^ is the umbral coun- 
terpart of a compound Poisson random variable with parameter 1. More in general, compound 
Poisson random variables of parameter t are represented by the auxiliary umbrae t./S.n^, with 
generating function f(t.f3.Ko,,z) = exp[t{f{Kc, z) — 1)] and we have t.a = t./S.K^. 

Let {Xt}t>o be a real-value Levy process, i.e. a process starting from and with stationary 
and independent increments. If we denote the moment generating function of Xt+s — Xg by 
<j){z,t), then (p{z,t) is infinitely divisible [12] and (j){z,t) = {(p{z,l)Y = (V'(-z))*, where tp{z) is 
the moment generating function of Xi. In particular, (l){z,t) = exp[t log ^l^{z)] = exp[tk{z)], where 
k{z) is the cumulant generating function of Xi such that A;(0) = 0. Comparing exp[tk{z)] with 
exp[t{f{Ka, z) — 1)], the correspondence between t.a = t./S.K^ and the Levy process {^j}j>o is 
immediate. We call t.a the Levy umbra associated to the umbra a. 

The following theorem explicitly states the connection between time-space harmonic polyno- 
mials with respect to Levy processes {Xt}t>o and the sequence of cumulants of Xi. 

Theorem 4.6. For all nonnegative integers k, the family of polynomials 

Qk{x, t) = E[{x - t.p.K^f] e (4.2) 
is time-space harmonic with respect to {t.a}t>o, where is the a-cumulant umbra. 

The umbra —t.a = t.{—l.a) is the Levy umbra associated to the umbra —l.a. Therefore, 
also the polynomials Qk{x,t) = E[{x + t.ji.K^)^] are time-space harmonic with respect to Levy 
processes umbrally represented by {—t.a}t>o- For simplicity, the following results refer to this last 
class of polynomials, but they can be stated also for the polynomials given in (j4.2p . 

Moments of a polynomial umbra as x + t.p.^ are generalizations of exponential complete Bell 
polynomials, see [Ij. 

Proposition 4.7 (Sheffer identity with respect to t). For the sequence of polynomials {Qk{x-,t)} 
in 14-'^ the following identity holds 

Qk{x,t + s) = Y,[ .]Pj{s)Qk-j{x,t), 
j=0 ^-^^ 

where Pj{s) = Qj{0,s) for all nonnegative integers j. 

Proof The result follows as Qk{x,t + s) = E {[x + t.p.K^ + s./S.k^]^) = Y.^-^o [{s.p.K^y] 
E[{x + t.l3.K^)''-^]. □ 
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Corollary 4.8. For the coefficients {qj''\t)} of the sequence of polynomials {Qk{x,t)} in 



we have 



i=l ^ ' 



for j = 1, . . . ,k, where {hi} is the sequence of cumulants of a. 

Proof. Consider t as an indeterminate and observe that ^qf^\t) = \{t.f5.Ka)^~^^^ . We have 

^'zf w - [[{t+x).p.^.r^ - {t.p.K.r^} ^ Q) g {^-'y\{t.p.K^f-^-\ (4.3) 

see [1] for the first equivalence. The result follows by taking the evaluation of both sides in 
(1431) . □ 

4.3 Special families of polynomials 

In this section we show that some classical families of time-space harmonic polynomials are mo- 
ments of the polynomial umbra a;— t. a, or can be recovered as a linear combination of its moments, 
for a suitable umbra a. 

Hermite polynomials. Standard Brownian motion is a special Levy process. From the 
Levy-Khintchine formula [12j, its generating function is (j){z., t) = exp(tzV2). Therefore, the 
umbral counterpart of a standard Brownian motion is the umbra where /(5, z) = 1 + j1. 

From Theorem l4.61 time-space harmonic polynomials with respect to a standard Brownian motion 
are Qfc(x,i) = E\{x - t.^.Sf\. 

Proposition 4.9. For all nonnegative integers k we have Qk{x,t) = Hj^\x), where Hjf \x) are 
the generalized Hermite polynomials with generating function X]^>o Hj^ — expjxz — s^}. 



(s ) 

Proof. We have Hjf (x) = E {x — l.p.{s6))'' , see [5]. The result follows by observing that 
-l.^.{Vi6) = -t.p.5. □ 

Poisson-Charlier polynomials. It is known that the Poisson-Charlier polynomials Ck{x + 
t, t) are time-space harmonic with respect to the compensated Poisson process [TT], with {Ck{x, t)} 
polynomials having generating function X^^>o Ck{x, t)j^ = e~*^(l+z)^'. We will recover this result, 

by proving that Ck{x + 1, t) is a linear combination of the polynomials Qkix, t) in p.2p when they 
are time-space harmonic with respect to the umbral counterpart of a compensated Poisson process. 

Let {Nt}t>o be a Poisson process of intensity 1 and Xt = Nt — t the compensated process 
|13j . From the Levy-Khintchine formula [12], its generating function is (j){z,t) = exp[t(e^ — z)]. 
The umbra 7 having generating function 7(7, z) = e^ — z is the disjoint difference of the unity 
umbra u and the singleton umbra x- Indeed recall that given two distinct umbrae 61 and 82, their 
disjoint difference is an auxiliary umbra, denoted by the symbol 61 — 82, with generating function 
f{6i,z) — f{52,z) + 1. Therefore, the umbral counterpart of a compensated Poisson process is 
t.f3.{u—x)- From Theorem 14.61 time-space harmonic polynomials with respect to a compensated 
Poisson process are Qk{x,t) = E[{x — t.(3.{u—x)}^]- 
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Proposition 4.10. We have Ck{x + t,t) = j)Qj(a;, i), where s{k,j) are the Stirling 

numbers of the first kind. 

Proof. By comparing generating functions, the sequence of polynomials {Ck{x, t)}k>o is umbrally 
represented by the polynomial umbra x.x — t. Therefore, we have Ck{x + t, t) = E[{{x + t).x — t)^]. 
We also have {x + t).x — t = {x.x — t)+ t.x = {{x.x — t)-/3 + t).x, due to the distributive property 
dTi]) and (3.x = u. Since E[{a.x)''] = E[{a)k], see [2], we have Ck{x + t,t) = E[{{x.x-t).(3 + t)k] = 
J2^j=osik,j)E[{{x.x-t).(3+ty]. The result follows by showing that x-t./3.(u-x) = {x.x-t).f3+t. 
Indeed we have x — t.f5.{u—x) = x — t.fi — t.(3.{—x) = {x — t.f3 — t./3.{—x))-X-f^ = {^-X — t — 
x)-x)-/5; where we have used again x./3 = u and the distributive property (j2.4p . The result 
follows since = X-(-l) and -t./3.{-x).X = -^•/3-X-(-l)-X = t-X- □ 

Levy-ShefTer polynomials. According to the definition given in [Ijj, a sequence of poly- 
nomials {Vk{x,t)}t>o is a Levy-Sheffer system if it is defined by the following generating function 

^k{x, t)^ = {g{z)Y exp{xu{z)}, (4.4) 

A:>0 

where g{z) and u{z) are analytic in a neighborhood of z = 0, u{0) = 0, ^(0) = 1, u'(0) ^ and 
l/g(r(z)) is an infinitely divisible moment generating function, with t{z) such that t{u{z)) = z. 
Schoutens in [14] states that the basic link between these polynomials and Levy processes is a 
martingale equality (cf. pag. 337 eq. (6)), which is equivalent to ask that these polynomials are 
time-space harmonic with respect to Levy processes. 

In the following, we show that Vk{x,t) is a linear combination of suitable time-space harmonic 
polynomials Qk{x,t) and therefore, they share the same property. 

To this aim we need to recall the notion of compositional inverse of an umbra a. Recall first that 
the auxiliary umbra a.f3.j is the composition umbra of a and 7, see [2]. The name recalls that its 
generating function is the composition of /(a, z) and /(7, z), that is /(a./3.7, z) = f{a, 7(7, z) — l). 
Moreover its moments are 

i 

). (4.5) 

k=l 

If E\a] 7^ 0, the compositional inverse of an umbra a is the auxiliary umbra a<~^> such that 
a.p.a^^^^ = Q^~^^./3.a = X) with x the singleton umbra. 

Theorem 4.11. We have Vk{x,t) = Ya=qE[{x + t.p.Ky]Bk,i{gi, . . . ,gk-i+i), where gj = E[-f^], 
for all nonnegative j and k is the cumulant umbra of a./S.'j^"^^ , with 7<"^> the compositional 
inverse of the umbra 7. 

Proof. Due to the form of the generating function in (j4.4p . the polynomials {Vkix,t)}k>o are 
umbrally represented by the polynomial umbra t.a + x./3.j, with f{a,z) = g{z) and f{'^.,z) = 
1 -|- u{z), and E['y\ = u'{0) ^ 0. Thanks to the distributive property (j2.4p we have x.f3.^ + t.a = 
(x -|- t.a./3.7^"^^)./3.7. The result follows by using (14. 5|) . □ 

Corollary 4.12. The Levy-Sheffer polynomials {Vk{x, t)}t>o are time-space harmonic with respect 
to Levy processes umbrally represented by {—t.a.(3.j'^~^^}t>o. 
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Time-space harmonic polynomials in terms of boolean and free cumulants. Let 

M{z) be the ordinary generating function of a random variable X, that is M{z) = 1 + ^j>x (liZ^, 
where Oj = We have M{z) = 1/[1 — B{z)], where B{z) = X^j>i hz^, and bi are the boolean 

cumulants of X. The noncrossing (or free) cumulants of X are the coefficients rj of the ordinary 
power series R{z) = 1 + Ylii>i''^i^^ such that M{z) = i?[zM(z)]. The umbral theory of boolean 
and free cumulants has been introduced in [6j. In particular, the a-boolean cumulant umbra 
r/„ has moments E[ri^^ = hi for all nonnegative integers i. This umbra is such that a = u./S.fj^, 
where = n! and £'[«"'] = nla^ for all nonnegative integers n. Thanks to Theorem 13.51 the 

polynomials Qk{x, t) = E[{x — t.u.(3.fja)''] are time-space harmonic polynomials with respect to the 
family {t.a}t>o with rj^ the a-boolean cumulant umbra. The a -free cumulant has moments 
= ilri, for all nonnegative integers i. This umbra allows us a different parametrization 
of Qk{x,t) = E[{x — t.a)'']. Indeed, if we denote by an the derivative umbra of a, such that 
f{ao,z) = 1 + zf{a,z), then we have a = .^a./3.(— LJI^)^"^^. Therefore, also the polynomials 
Qk{x,t) = E[{x + 1..^q)./3.(— I.jIq)^"^^)'^] are time-space harmonic polynomials with respect 
to the family {t.a}t>o- 
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